A0-fti33  OH3 

Special  Report  83-22 

July  1983 


US  Army  Corps 
of  Engineers 

Cold  Regions  Research  & 
Engineering  Laboratory 


Estimating  transient  heat  flows  and 
measuring  surface  temperatures  of 
a  built-up  roof 


Charles  J.  Korhonen 


>- 

a- 

o 

CJ> 


UJ 


DTIC 


Approved  for  public  release;  distribution  unlimited 


Unclassified 


SECURITY  CLASSIFICATION  of  This  PACE  (Whtn  Data  Entaram 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

Special  Report  83-22 

WKnSEBGn 

3.  RECIPIENT'S  CATALOG  NUMBER 

«.  TITLE  One!  Subtitle) 

ESTIMATING  TRANSIENT  HEAT  FLOWS  AND 

MEASURING  SURFACE  TEMPERATURES 

OF  A  BUILT-UP  ROOF 

S.  TYPE  OF  REPORT  &  PERIOD  COVERED 

S-  PERFORMING  ORG.  REPORT  NUMBER 

T.  AUTHOR!  a) 

Charles  J.  Korhonen 

a.  CONTRACT  OR  GRANT  NUMBER!.) 

»-  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

U.S.  Army  Cold  Regions  Research  and 

Engineering  Laboratory 

Hanover,  New  Hampshire  03755 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

IS.  NUMBER  OF  PAGES 

21 

14.  MONITORING  AGENCY  NAME  A  ADDRESS!!/  dlllarant  from  Controlling  Ottlca) 

IS.  SECURITY  CLASS,  (ot  thla  r.porf) 

Unclassified 

1S«.  DtCLASSlFICATION/ DOWN  GRADING 

schedule 

t«.  DISTRIBUTION  STATEMENT  (of  Oil  a  Report) 


Approved  for  public  release;  distribution  unlimited. 


IT.  DISTRIBUTION  STATEMENT  (of  U >•  abatract  antarad  |n  Bloc*  29,  It  dlllarant  turn,  Report) 


IS.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  rmveree  tide  If  neceeeery  end  Identify  by  block  number) 

Heat  Roofs 

Heat  flow  Thermovision 

Infrared  equipment 
Roof  temperature 

'  20.  ABSTRACT  (Cmntbm*  an  nnm  mid*  It  nag.  ..ary  arW  Idaotlty  by  block  ni mibar) 

transient  heat  flow  through  a  multilayered  building  component  can  be  estimated 
using  the  transfer  function  method  presented  in  the  ASHRAE  (1977)  Handbook  of 
Fundamentals .  Sol-air  temperature  is  one  parameter  recommended  for  use  in  this 
method,  but  surface  temperatures  were  shown  to  be  a  reasonable  substitute.  Al¬ 
though  the  magnitude  of  the  heat  flow  as  calculated  with  the  transfer  function 
appears  to  be  reasonable,  more  testing  should  be  carried  out  to  determine  its 
accuracy.  An  infrared  camera  can  measure  roof  surface  temperatures  fairly  ac- 
curately;  the  most  accurate  measurements  were  made  at  night,. _ 

DO  iISTt,  U73 


EDITION  OF  •  NOV  as  IS  OBSOLETE 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PACE  Data  Entarad) 


PREFACE 

This  report  was  prepared  by  Charles  J.  Korhonen,  Research  Civil 
Engineer,  Civil  Engineering  Research  Branch,  Experimental  Engineering 
Division,  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory,  while 
he  was  at  the  University  of  Alaska,  Fairbanks,  from  September  1980  to  June 
1981.  This  report  was  submitted  to  the  author's  degree  committee  as  part 
of  the  requirement  for  completion  of  an  advanced  degree  in  Arctic  Engineer¬ 
ing.  The  degree  committee  consisted  of  John  L.  Burdick,  Head,  Department 
of  Civil  Engineering;  Elbert  F.  Rice,  Professor  of  Civil  Engineering;  and 
John  P.  Zarling,  Head,  Department  of  Mechanical  Engineering. 

Technical  review  was  provided  by  S.  Flanders  and  G.  Phetteplace  of 
CRREL. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  pro¬ 
motional  purposes.  Citatiou  of  brand  names  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  of  such  commercial  products. 
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ESTIMATING  TRANSIENT  HEAT  FLOWS  AND  MEASURING  SURFACE 
TEMPERATURES  OF  A  BUILT-UP  ROOF 
by 

Charles  J.  Korhonen 


INTRODUCTION 

To  properly  size  heating  and  air-conditioning  systems  for  buildings, 
it  is  important  to  know  the  hourly  heat  flows  through  the  building's 
envelope.  Because  building  elements  not  only  resist  the  flow  of  heat  but 
store  it  as  well,  the  maximum  heating  or  cooling  demand  in  the  building  is 
not  usually  In  phase  with  the  cyclic  changes  in  ambient  temperatures.  A 
sudden  increase  or  decrease  in  temperatures  across  a  building  element  will 
not  cause  a  corresponding  change  In  heat  flow  through  the  element.  The 
ASHRAE  (1977)  Handbook  of  Fundamentals  presents  a  calculation  procedure 
that  an  engineer  can  use  to  estimate  instantaneous  rates  of  heat  loss-gain 
by  knowing  certain  facts  about  the  building  envelope  and  the  weather. 

My  first  purpose  was  to  verify  this  technique  by  measuring  actual  heat 
flows  through  a  model  built-up  roof  with  a  heat  flow  sensor.  My  second 
purpose  was  to  measure  the  surface  temperatures  on  that  model  roof  with  an 
infrared  scanner  and  a  thermocouple.  Since  infrared  scanners  are  becoming 
popular  as  a  tool  for  diagnosing  building  temperatures  and  heat  flow 
problems,  a  comparison  was  made  between  the  Infrared  scanner  and  thermo¬ 
couple  results. 

EXPERIMENT  SETUP 

A  A-  by  4-ft  model  built-up  roof  was  designed  and  built  to  serve  as 
the  building  component  on  which  to  evaluate  the  ASHRAE  (1977)  Handbook  of 
Fundamentals  computation  technique  and  on  which  to  measure  surface  tempera¬ 
tures,  TTie  model  roof  was  placed  on  an  electrically  heated  box  that  was 
set  f*n  the  roof  of  the  Duckerlng  Building  on  the  Fairbanks  campus  of  the 
University  of  Alaska.  Each  layer  of  the  model  roof  was  carefully  dimen¬ 
sioned  and  assembled  to  provide  a  building  component  of  known  cross 
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Figure  1.  Cross  section  of 
built-up  roof. 


Figure  2.  Insulation  pattern 
in  built-up  roof. 


section.  Thermocouples ,  placed  between  each  layer,  and  a  heat  flow  sensor, 
taped  to  the  underside  of  the  roof,  were  read  hourly  from  3  a.m.  to  12  mid¬ 
night  on  24  April  1981  to  furnish  data  for  evaluating  the  Handbook  of 
Fundamentals  technique  (Appendix  A).  A  thermocouple  was  embedded  just 
below  a  film  of  roof  cement  to  measure  surface  temperatures. 

Figure  1  shows  the  arrangement  of  the  heat  flow  sensor  and  copper- 
constantan  thermocouples.  The  insulation  in  the  roof  consisted  of  four 
1-in.  layers  of  perlite  boardstock,  which  is  a  blend  of  perlite  particles 
and  reprocessed  cellulose  fibers.  A  3/4-in.  plywood  sheet  served  as  the 
deck.  The  built-up  membrane  consisted  of  alternate  layers  of  asbestos 
cement  that  was  applied  cold  and  three  plies  of  no.  If  roofing  felt.  No 
gravel  was  applied  to  the  roof  surface.  (As  a  note  of  interest,  although 
care  was  taken  to  exclude  all  voids  during  the  construction  of  the  built-up 
membrane,  a  few  small  blisters  yet  developed  in  the  membrane.  It  is  easy 
to  understand  how  blisters  can  be  built  into  actual  roofs.) 

Figure  2  shows  the  layout  of  the  2-ft-wide  insulation  boards  in  the 
4-ft-wide  roof.  This  pattern  created  a  2-  x  2-ft  test  area  for  measuring 
heat  flows  that  was  free  of  insulation  seams. 

Five  60-W  light  bulbs  connected  to  a  thermostat  provided  heat  whenever 
the  box  temperature  dropped  below  75°F.  The  inside  air  was  continually 
mixed  by  a  small  electric  fan.  The  heat  flow  sensor  was  shielded  from  the 
light  bulbs  by  a  sheet  of  plywood. 
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RESULTS 

Transfer  function 


The  transfer  function  method  is  presented  in  Handbook  of  Fundmentals 
(ASHRAE  1977)  for  calculating  the  heat  conducted  through  multilayered 
building  elements  exposed  to  transient  thermal  conditions.  The  Handbook 
recommends  that  "sol-air"  temperatures  he  used  in  carrying  out  the  calcula¬ 
tions,  but  as  I  will  show,  the  use  of  roof  surface  temperatures  yields 
similar  results.  Thus,  it  might  be  possible  to  use  an  infrared  scanner  to 
remotely  gather  these  data. 

The  sol-air  temperature  can  be  represented  by  the  following  equation: 

T  =  T  +  «  i  /h  -  eAR/h  (1) 

e  o  t  o  o 

where  Te  =  sol-air  temperature  (°F) 

T0  =  outdoor  air  temperature  (°F) 

*  =  solar  absorptance 
It  «  incident  solar  radiaton  (Btu/hr*ft2) 
hg  =*  air  film  heat  transfer  coefficient  (Btu/hr*ft  *°F) 
e  *  surface  emissivity 

AR  -  difference  between  the  longwave  radiation  incident  on  the  sur¬ 
face  from  the  sky  and  surroundings,  and  the  radiation  emitted 
by  a  black  body  at  outdoor  air  temperature  (Btu/hr *ft2) . 

An  estimate  of  the  solar  radiation  incident  on  the  test  roof,  Tt ,  is 
given  for  latitude  64  degrees  north  in  Table  25  of  the  Handbook  of  Funda¬ 
mentals  (ASHRAE  1977,  Chapter  26).  The  listed  values  cannot  be  used 
directly  because  they  represent  the  solar  heat  gain  through  double-strength 
sheet  glass.  Consequently  they  must  first  be  multiplied  by  1.15  (ASHRAE 
1977)  to  correct  for  the  solar  energy  excluded  by  the  glass.  Since 
Fairbanks  was  on  daylight  savings  time  during  this  experiment,  the 
tabulated  data  were  shifted  l  hour.  The  corrected  and  shifted  values  of 
It  from  the  Handbook  are  listed  in  Table  1,  These  values  increase  from  a 
very  small  thermal  load  in  the  morning  to  a  maximum  at  1  p.m.  (solar  noon) 
and  then  decrease  in  like  manner  in  the  afternoon.  The  day  on  which  this 
experiment  was  run  was  nearly  clear  so  no  corrections  were  made  for  cloud 
cover . 
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Table  2.  Sol-air  and  test  roof  surface 
temperatures. 


Sol-air  Roof  surface 

Time  (°F) _ (°F) 


1  a.m. 

31.2 

31 .8 

2 

29.4 

27.2 

Table  1.  Solar 

radiation 

3 

27.6 

22.6 

incident  on  test 

roof  (clear 

4 

29.2 

24.8 

day). 

5 

22.0 

29.0 

6 

31.9 

36.2 

Local  time 

I 

7 

44.1 

59.0 

(hr) 

(Btu?hr .ft2) 

8 

55.2 

66.8 

9 

76.5 

109.8 

6  a.m. 

2.3 

1.0 

89.4 

111.4 

7 

24.2 

11 

100.9 

121.2 

8 

58.7 

Noon 

112.1 

121.2 

9 

9V.8 

1  p.m. 

111.0 

103.4 

2 

110.9 

109.4 

11 

161.0 

3 

103.7 

82.4 

Noon 

178.3 

4 

92.4 

73.0 

1  p.m. 

184.0 

5 

81.3 

66.4 

2 

178.3 

6 

69.4 

60.8 

3 

161,0 

7 

52.7 

46.6 

4 

133.4 

8 

42.1 

39.0 

5 

97.8 

9 

39.6 

35.6 

6 

58.7 

10 

39.0 

36.0 

7 

24.3 

11 

38.2 

37.8 

8 

2.3 

12 

33.0 

36.4 

The  Handbook  indicates  that  a  flat  roof  surface  loses  more  longwave 
radiation  than  it  gains  from  the  surroundings.  As  a  consequence,  eAR/h0 
equals  -7°?. 

By  measuring  air  temperatures  and  selecting  e  and  hQ  values,  the 
engineer  can  use  eq  1  to  calculate  sol-air  temperatures.  For  example,  the 
aol-air  temperature  at  10  a.m.  on  24  April  with  e  ■  0.93  and  h0  =  3  is 

0  °3 

T  -  55*  +  (133.4)  -  7 

e  j 

T  =*  89.4eF. 
e 

The  roof  surface  temperature  (Table  2)  at  that  same  hour  was  lll.4°F. 
Although  this  is  not  good  agreement,  the  sol-air  and  roof  temperatures  were 


*  The  air  temperatures  are  taken  from  Appendix  A. 


Table  3.  Heat  flows  (3tu/hr)  through  test  roof  calculated  with  eq  2  using 
sol-air  temperatures. 


Time 

(hr) 

First 

iteration 

Se  cond 
iteration 

Third 

iteration 

Fourth 

iteration 

1 

-0.9696940 

-2.4442700 

-2.4442710 

-2.4442710 

2 

-1.8340240 

-2.7049590 

-2.7049600 

-2.7049600 

3 

-2.4379300 

-2.9298360 

-2.9298360 

-2.9298360 

4 

-2.8442920 

-3.1179190 

-3.1179190 

-3.1179190 

5 

-3.1046800 

-3.2560730 

-3.2560730 

-3.2560730 

6 

-3.2949180 

-3.3785220 

-3.3785220 

-3.3785220 

7 

-3.4269170 

-3.4730540 

-3.4730540 

-3.4730540 

8 

-3.2887370 

-3.3141920 

-3.3141920 

-3.3141920 

9 

-2.8673470 

-2.8813900 

-2.8813900 

-2.8813900 

10 

-2.1868970 

-2.1946430 

-2.1946430 

-2.1946430 

11 

-1.2659210 

-1.2701940 

-1.2701940 

-1.2701940 

12 

-0.2885270 

-0.2908844 

-0.2908843 

-0.2908843 

13 

0.6516642 

0.6503637 

0.6503637 

0.6503637 

14 

1.4486700 

1.4479530 

1. 4*79530 

1.4479530 

15 

1.9694940 

1.9690980 

1.9690980 

1.9690980 

16 

2.2109110 

2.2106920 

2.2106920 

2.2106920 

17 

2.1383280 

2.1382070 

2.1382070 

2.1382070 

18 

1.7770200 

1.7769540 

1.7769540 

1.7769540 

19 

1.2124580 

1.2124210 

1.2124210 

1.2124210 

20 

0.4780363 

0.4780160 

0.4780160 

0.4780160 

21 

-0.3839333 

-0.3839445 

-0.3839445 

-0.3839445 

22 

-1.1876650 

-1.1876710 

-1.1876710 

-1 . 187b7 10 

23 

-1.7802520 

-1.7802560 

-1.7802560 

-1.7802560 

24 

-2.1646280 

-2.1646300 

-2.1646300 

-2.1646300 

closer  during  the  evening  hours.  Since  the  sol-air  temperatures  are  based 
on  estimated  incoming  solar  radiation,  and  roof  surface  temperatures 
reflect  actual  conditions.  I  feel  that  the  thermocouple  readings  will  yield 
the  best  results  in  the  transfer  function  calculations. 

Equation  2,  which  is  recommended  by  the  Handbook  for  estimating  hourly 
heat  flows  through  a  building  element,  was  used  to  calculate  the  heat  flow 
values  in  Tables  3  and  4.  The  sol-air  and  surface  temperature  data  in 
Table  2  simulated  the  outside  air  temperatures.  Inside  air  temperatures 
were  held  at  75° F,  while  the  indoor  and  outdoor  air  film  coefficients  were 
the  same  as  those  used  in  eq  1. 


e,r 

where 


q  -  nA 

Alt  b  (t  -  nt)  -  Z  d  - )  -  t  1  c  ] 

n-0  n=*l  n-0 

A  *  indoor  surface  area  (ft2) 
qe  T  ■  heat  flow  through  indoor  surf  ice  (Btu/hr) 


(2) 
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Table  4.  Heat  flows  (Btu/hr)  through  test  roof  calculated  with  eq  2  using 
actual  roof  temperatures. 


Time 

(hr) 

First 

iteration 

Second 

iteration 

Third 

iteration 

Fourth 

iteration 

1 

-0.9573820 

-2.6242310 

-2.6242320 

-2.6242320 

2 

-1. 7733UO 

-2.7543960 

-2.7543960 

-2.7543960 

3 

-2.3911340 

-2.9446140 

-2.9446140 

-2.9446140 

4 

-2.8852500 

-3.1930050 

-3.1930050 

-3.1930050 

5 

-3.2545280 

-3.4247800 

-3.4247800 

-3.4247800 

6 

-3.4288710 

-3.5228840 

-3.5228840 

-3.5228840 

7 

-3.3901210 

-3.4420010 

-3.4420010 

-3.4420010 

8 

-3.0390340 

-3.0676570 

-3.0676570 

-3.0676570 

9 

-2.3527220 

-2.3685120 

-2.3685130 

-2.3685130 

10 

-1.3356600 

-1.3443700 

-1.3443700 

-1.3443700 

11 

0.0230660 

0.0182607 

0.0182607 

0.0182607 

12 

1.1916790 

1.1890280 

1.1890280 

1.1890280 

13 

2.0744340 

2.0729720 

2.0729720 

2.0729720 

14 

2.5047830 

2.5039760 

2.5039760 

2.5039760 

13 

2.5186090 

2.5181640 

2.5181640 

2.5181640 

16 

2.2870190 

2.2867740 

2.2867740 

2.2867740 

17 

1.6627210 

1.6625860 

1.6625860 

1.6628860 

18 

0.9325699 

0.9324952 

0.9324952 

0.9324952 

19 

0.2697038 

0.2696626 

0.2696626 

0.2696626 

20 

-0.3697135 

-0.3697362 

-0.3697362 

-0.3697362 

21 

-1.0655870 

-1.0655990 

-1.C655990 

-1.0655990 

22 

-1.7079740 

-1.7079810 

-1.7079810 

-1.7079810 

23 

-2.1906630 

-2.1906670 

-2.1906670 

-2.1906670 

24 

-2.4769670 

-2.4769690 

-2.4769690 

-2.4769690 

t  *  time 

A  *  time  Interval 

n  **  number  of  transfer  function  coefficients 
te,x  ~  ■  sol-air  or  roof  temperature  at  time  t  _  nA  (°F) 

trc  “  indoor  temperature  (°F) 
bn,dn,cn  *  transfer  function  coefficients. 

The  transfer  function  coefficients  d>n»^n'cn^  are  Hsted  for 
typical  wall  and  roof  constructions  in  the  Handbook.  Alternatively,  the 
engineer  can  use  Mitalas  and  Arseneault's  (1972)  computer  program  to 
calculate  transfer  function  coefficients,  as  I  did  for  this  report.  The 
University  of  Alaska  has  this  program  on  file  in  its  computer.  The 
computer  calculated  the  transfer  function  coefficients  after  I  supplied  it 
with  an  outdoor  film  coefficient  h0  ■  3.0  Btu/hr»ftz»°F  and  an  indoor 
fjlm  coefficient  -  1.46  Btu/hr»ft2»°F,  and  the  density,  conductivity, 
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specific  heat,  thickness,  emissivity  and  absorptivity  of  the  roof. 

The  resulting  tranfer  function  coefficients  for  the  test  roof  are: 

b0  -  0.00018  d0  -  1.00000 

bj  -  0.00795  =  -0.79820 

b2  -  0.01409  I  Cn  =  0.02526  d2  -  0.14593 

b3  -  0.00297  n-0  d3  =■  -0.00552 

bu  -  0.00007  d4  -  0.00003 

Inspection  of  eq  2  reveals  that  the  calculation  of  heat  flow  for  any 
given  hour  requires  temperatures  at  that  hour  and  preceding  hours,  plus 
heat  flows  at  preceding  hours.  Since,  initially,  heat  flows  are  not  known, 
any  value  of  heat  flow  can  be  used  as  a  first  guess.  Through  an  iterative 
process  the  actual  calculated  heat  flow  will  evolve. 

Equation  2  is  composed  of  three  sections.  The  first  section  sums  the 
product  of  each  b  coefficient  and  simulated  outdoor  condition.  The  present 
outdoor  temperature,  either  sol-air  or  roof  surface,  is  multiplied  by  the 
b0  coefficient.  Next  the  outdoor  temperature  from  the  previous  hour  is 
multiplied  by  the  b3  coefficient.  This  continues  until  all  the  b  coef¬ 
ficients  are  used  and  then  the  products  are  summed  to  complete  the  first 
section  of  eq  2. 

The  second  section  is  the  summation  of  the  d  coefficients  multiplied 
by  the  previous  values  of  heat  flow.  Equation  2  does  not  use  d0  in  its 
calculations.  The  dj  coefficient  is  multiplied  by  the  heat  flow  calculated 
for  the  previous  hour,  then  d2  is  multiplied  by  the  heat  flow  2  hours  back 
in  time.  This  continues  until  all  the  d  coefficients  (except  dQ)  are 
used.  The  products  are  then  summed  to  complete  the  second  section  of  eq 
2.  The  third  section  is  calculated  only  once  since  there  is  only  one  c 
coefficient  and  the  inside  temperature  is  constant. 

Once  all  three  sections  are  calculated  they  are  combined  to  give  the 
values  of  heat  flows  for  each  hour  of  the  day.  This  process  is  repeated 
until  the  latest  calculated  heat  flow  closely  agrees  with  the  previous 
iteration.  This  convergence  of  values  is  shown  in  Tables  3  and  4.  In 
this  case  either  the  third  or  fourth  iteration  yields  the  desired  estimated 
heat  flow.  Appendix  B  lists  the  BASIC  computer  program  written  to  perform 
the  calculations. 
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Heat  flow  and  temperature  measurements 


Copper- conatantan  thermocouples  were  used  for  this  experiment.  A 

Fluke  Model  2175A  Digital  Thermometer  was  used  to  measure  the  thermocouple 

readings.  It  maintained  its  own  internal  reference  and  automatically  con¬ 
verted  the  millivolt  signal  from  the  thermocouples  into  a  temperature  In 
degrees  Fahrenheit, 

The  heat  flow  sensor  used  in  this  experiment  consisted  of  a  number  of 
thermocouples  embedded  in  a  material  of  known  thermal  conductivity.  Half 
of  the  thermocouples  were  embedded  near  the  top  surface  and  half  near  the 
bottom  surface  of  the  sensor  to  measure  the  temperature  change  across  the 
sensor.  This  produced  a  millivolt  reading  proportional  to  heat  flow.  This 
heat  flow  sensor  had  a  calibration  constant  of  93.1  mV,  equal  to  500 
Btu/hr»ft2.  A  Keithley  Instruments  Model  155  Microvoltraeter  was  used  to 
obtain  heat  flow  sensor  readings.  Data  obtained  from  the  heat  flow  and 
temperature  sensors  are  presented  in  Appendix  A. 


Infrared  temperature  measurements 

I  used  an  Infrared  camera  to  determine  how  accurately  surface  tempera¬ 
tures  could  be  measured.  An  AGA  Thermovision  750,  which  measured  radiant 
energy  within  a  wavelength  band  of  2  to  5.6  pm,  was  used  for  these  measure¬ 
ments. 

The  energy  emitted  from  the  surface  of  a  body  is  a  function  of  the 
object' 8  temperature  and  emissivity.  Emissivity  relates  to  the  ease  with 
which  cn  object  emits  infrared  radiation.  Not  all  surfaces  emit  the  same 
amount  of  energy  when  they  are  at  identical  temperatures.  A  surface  that 
emits  the  maximum  amount  of  energy  has  an  emissivity  of  1.  All  objects 
In  nature  possess  emlssivltles  less  than  i.  Therefore,  an  accurate 
measurement  of  surface  temperature  requires  a  knowledge  of  the  longwave 
emissivity  of  the  object  to  correct  for  varying  radiation  potentials. 

Equation  3  accounts  for  both  radiation  emitted  and  reflected  from  the 
surface  of  an  object.  (This  emissivity  equation  is  presented  in  the 
owner's  manual  [AGA  Infrared  Systems  1973]  for  the  infrared  camera.) 
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where  Ai 


or  =*  image  isotherm  difference  (i0  -  ir)  read  from  display 
unit 

eQ  =  object  emissivity  (roof) 
er  =  reference  emissivity  (timber) 

Ir  ■  absolute  isotherm  read  from  calibration  chart  (Fig.  3) 
knowing  reference  temperatures  Tr 
la  =  absolute  isotherm  read  from  calibration  chart  (Fig.  3) 
knowing  ambient  temperature  Ta. 

By  knowi  t  the  temperature  and  emissivity  of  a  reference  object  and  the  air 
temperat  2  and  the  emissivity  of  another  object  we  can  use  eq  3  with 
Figure  o  calculate  the  unknown  temperature  of  the  second  object. 

I  measured  radiation  with  an  old  wooden  timber  of  an  assumed  emis¬ 
sivity  of  0.9  in  the  same  field  of  view  as  the  model  roof.  One  thermo¬ 
couple  monitored  the  timber's  surface  temperature  while  a  second  measured 
the  air  temperature.  (The  air  temperature  thermocouple  was  shaded  from  the 
sun.)  The  roof  was  assumed  to  have  a  high  emissivity  of  0.93  since  it  was 
fresh,  clean  asphalt. 

During  one  measurement  cycle  the  timber  and  air  temperatures  were 
44.6°F  and  46.9°F  respectively.  The  resulting  absolute  isotherm  levels 
I_  and  I3  as  read  from  a  calibration  chart  (Fig.  3)  were  12.8  and  13.4 

l  a 

respectively.  The  isotherm  markers  of  the  timber  and  the  roof,  read  from 
the  camera  display  unit,  were  0.5  and  0,4  respectively.  The  image  isotherm 
levels  ir  and  i0  were  then  obtained  by  multiplying  the  isotherm  markers 


T«mp«rotur*  (*f)  Figure  3.  Calibiation  chart,  f/1.8 

denotes  aperture  setting  of  camera. 
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Table  5.  Surface  temperature  of  test  roof  as  measured  by  a 
thermocouple  and  an  Infrared  camera  on  24  April  1981. 


Measurement 

number 

Thermocouple 
no.  7 
(°F) 

I.R. 

camera 
(  ®F) 

Time 

Conditions 

1 

112.2 

111.2 

10:30  a.m. 

Sunny 

2 

113.0 

106.7 

12:50  p.ra. 

Partly  cloudy 

3 

91.4 

83.8 

1:15  p .m. 

Cloudy 

4 

102.6 

95.9 

1:20  p.ra. 

Sunny 

5 

105.6 

98.6 

1:22  p.ra. 

Sunny 

6 

108.2 

102.4 

1:25  p.ra. 

Sunny 

7 

37.8 

35.2 

9:30  p.m. 

Clear 

8 

37.6 

36.5 

9:40  p.ra. 

Clear 

9 

37.4 

35.6 

9:45  p.ra. 

Clear 

37.6 

36.5 

9:50  p.ra. 

Clear 

11 

37.9 

36.9 

11:30  p.m. 

Clear 

by  the  range  setting  of  the  camera,  which  was  20  during  this  measurement 
cycle;  ir  and  i0  became  10  and  8  respectively.  IQ  was  then  deter¬ 
mined  by  inserting  the  above  values  into  eq  3; 


1 

o 


8-10  0.90 

0.93  0.93 


(12.8)  +  (1  -  £||)  13.4 


I  -  10.7. 
o 


The  roof  surface  temperature  taken  from  Figure  3  was  36.9°F.  Thermocouple 
7  (Fig.  1)  measured  the  roof  surface  temperature  to  be  37.9°F. 

I  repeated  this  procedure  a  number  of  times;  the  results  are  shown  in 
Table  5.  The  greatest  discrepancy  between  the  roof  temperatures  measured 
with  the  thermocouple  and  those  measured  with  the  infrared  camera  occurred 
during  the  daytime.  Measurement  3  in  Table  5  was  7.6°F  below  that  of  its 
corresponding  thermocouple  measurement.  The  evening  radiation  measurements 
were  all  within  3°F  of  the  thermocouple  temperatures.  The  best  correlation 
occurred  later  at  night.  In  all  cases  the  radiation  measurements  felL 
somewhat  below  the  thermocouple  measurements.  Although  Table  5  does  not 
show  complete  agreement  between  the  two  measurement  techniques,  it  does 
show  that  surface  temperatures  can  be  approximated  with  an  infrared  camera, 
especially  during  the  night. 


Uncertainty  analysis 

An  area  of  uncertainty  in  determining  unknown  temperatures  from  radia 
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tion  measurement  lies  in  how  accurate  the  values  substituted  into  eq  3 
are.  Holman  (1971)  describes  a  method  for  specifying  the  uncertainty  in 
the  results  of  experimental  data.  His  analysis  will  be  followed  to  deter¬ 
mine  the  uncertainty  in  the  results  of  eq  3. 

The  result,  I0,  of  eq  3  is  a  function  of  the  independent  variables 
Aior »  ^o*  er »  la*  thus 


I 

o 


(4) 


If  Wr  represents  the  uncertainty  of  IQ  then  Wj ,  W2,  W3,  W4  and  W$ 
represent  the  uncertainty  or  range  of  possible  values  of  the  independent 
variables.  By  specifying  this  range  of  possible  values  for  each  in¬ 
dependent  variable,  we  can  determine  the  uncertainty  in  I0  as  shown 
below: 


W 

r 


E(|s—  “‘>2  +  «2>2  -  fe2  “s)2 


or 


•3e 


3c 


31  31 

+  (tt  “«12  +  <sr  “=)21  •  (5> 

r  a 

The  independent  variables  have  the  following  ranges:  e0  “  0.93  ±  0.02, 
er  -  0.90  ±  0.02,  ea  *  1.0  t  0.02,  lr  *  12.8  ±  0.20,  Ia  -  13.4  ± 

0.20.  Although  Ir  and  Ia  are  functions  of  temperature,  the  chart  from 
which  Ir  and  Ia  are  taken  controls  the  uncertainty  in  these  values.  I 
feel  that  Ir  and  Ia  can  only  be  read  within  ±  0.2  of  their  actual 
values  as  shown  above.  Aior  is  a  function  of  how  well  the  isotherm 
marker  can  be  read  from  the  display  unit  of  the  camera.  Due  to  the  flicker 
of  the  cathode-ray  tube,  I  do  not  feel  that  I  could  consistently  read  the 
isotherm  marker  closer  than  ±  0.02, 

Using  eq  5  along  with  the  stated  ranges  of  the  independent  variable 

gives 

31  /All  =  1/e  -  1/0.93  =  1.08 

o  or  o  - 


31  /  3e 
o  o 


=  -Ai 


e  I  /e 
r  r 


2 

o 


+  e  I  /e 
r  a 


2 

o 


=  -  -2/0. 932  -  0. 30(12. 8)/0.932  +  0.90( 13.4)/0.932 
-  2.94 
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Figure  4.  Thermocouple  temperatures. 


31  /3e  -  I  U  -  I  / e 
or  r  o  a  o 


*  12.8/0.93  -  13.4/0.93 
»  -0.65 


31  / 3 I  -  e  /e  -  0.90/0.93  -  0.97 
or  r  o  - 


31/31  *  1  -  e  /e  -  1  -  0.90/0.93  =*  0.03 

o  a  r  o  - — 

The  uncertainty  in  measuring  I0  (eq  3)  becomes 

Wr  -  [(1.08  x  0.02 )2  +  (2.94  x  0.02)2  +  (-0.65  x  0.02)2 

+  (0.97  x  0.2  )2  +  (0.03  x  0.2  )2]1/2 

W  -  0.20 

r 

Thus  the  Ig  reading  of  10.7  from  the  Infrared  Temperature  Measurement 
section  obtained  with  the  infrared  camera  is  accurate  to  ±  0.20  unit. 


Evaluation  of  transfer  function 

The  transient  nature  of  the  heat  flow  through  the  test  roof  is  shown 
by  the  graph  of  thermocouple  temperatures  in  Figure  4  (values  taken  from 
Appendix  A).  The  effect  of  a  maximum  surface  temperature  occurring  at  11 

12 


Table  6. 

Heat  flow  through  test  roof 

(Btu/hr) 

A  negative 

sign  indicates  heat  being 

lost 

to  the  outdoors. 

Transfer  function 

(based  on  roof 

Heat  flow 

Time 

temperatures) 

sensor 

1  a.  ra. 

2 
3 


a.m.  Is  not  sensed  by  the  underside  of  the  roof  until  about  2  p.m.  This 
should  come  as  no  surprise  when  the  heat  storage  effects  of  the  insulation 
are  considered.  We  also  see  a  maximum  heat  flow  at  the  inside  roof  surface 
at  2  p.m. 

Does  the  transfer  function  predict  the  same  thermal  lags  shown  in 
Figure  4?  Table  6  compares  heat  flow  through  the  roof  as  computed  with  the 
transfer  function  to  that  measured  with  the  heat  flow  sensor,  (Roof  sur¬ 
face  temperatures,  which  are  actual  values,  were  used  to  derive  the  trans¬ 
fer  function  values  as  opposed  to  using  the  estimated  sol-air  tempera¬ 
tures.)  The  heat  flow  sensor  data  were  taken  from  Appendix  A.  Generally, 
the  transfer  function  predicts  the  same  thermal  lag  as  shown  in  Figure  4. 
Also,  the  periods  of  heat  loss  and  gain  correspond  for  the  results  from  the 
transfer  function  and  the  heat  flow  sensor.  Both  measurement  techniques 
show  heat  being  lost  from  the  roof  until  11  a.m.,  then  the  roof  gaining 


heat  until  6  or  7  p.w. ,  whereafter  it  again  began  to  lose  heat.  It  appears 
that  the  transfer  function  is  able  to  model  the  time  and  direction  of  heat 
flow. 

However,  the  magnitude  of  heat  flow  in  Table  6  differs  considerably 
between  the  transfer  function  results  and  the  heat  flow  sensor  results. 
Because  of  the  short  measurement  period  of  this  experiment,  it  is  difficult 
to  assess  this  discrepancy  other  than  to  explain  some  of  the  potential 
sources  of  error. 

In  reading  the  microvoltmeter  it  was  often  necessary  to  taka  an  aver¬ 
age  of  several  readings  to  get  one  hourly  reading  because  the  meter  needle 
often  swung  bade  and  forth.  At  first  I  thought  the  cycling  on  and  off  of 
the  lights  to  be  the  probable  cause  of  the  apparent  jumps  in  the  heat  flow 
readings,  but  careful  control  of  their  cycling  did  nothing  to  dampen  the 
swing  of  the  needle.  Even  turning  the  fan  off  did  not  help.  Consequently, 
I  used  a  system  of  averaging  several  readings.  Perhaps  some  other  sort  of 
recorder,  such  as  a  stripchart  recorder,  would  have  yielded  better  results. 

SUMMARY  AND  CONCLUSIONS 

I  used  a  computer  method  that  was  modified  to  employ  roof  surface 
temperatures  and  heat  storage  capacities  of  materials  to  estimate  the 

hour-by-hour  transient  heat  flow  through  a  known  cross  section  of  built-up 
roof.  Thermocouples  and  a  heat  flow  sensor  provided  a  check.  Also,  I  com¬ 
pared  roof  surface  temperatures  measured  with  an  infrared  camera  to  roof 
surface  temperatures  measured  with  a  thermocouple. 

The  computer  method  predicted  the  time  and  the  direction  of  heat  flow 
through  the  built-up  roof  but  the  magnitude  of  those  heat  flows  did  not 
correspond  well  to  heat  flows  measured  by  the  heat  flow  sensor.  This  could 
be  due  to  the  insensitivity  of  the  heat  flow  sensor  and  to  the  instability 
of  the  microvoltmeter  needle  when  readings  were  taken.  The  accuracy  of  the 
transfer  function  should  be  further  checked  by  use  of  better  equipment  and 
longer  measurement  cycles. 

Surface  temperature  measurements  taken  with  the  infrared  camera  cor¬ 
responded  fairly  well  with  the  roof  surface  temperatures  recorded  by  a 
thermocouple.  I  found  the  best  results  at  night  when  there  was  no  solar 
influence. 
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APPENDIX  A:  THERMOCOUPLE  AND  HEAT  FLOW  READINGS 

Thermocouple  readings  in  °F;  heat  flow  readings  in  mV;  negative  sign 
Indicates  heat  flowing  from  inside  surface  of  roof  up  into  insulation 


:l  PRINT*********************  INFORMATION  ********************** * 

2  PRINT 

3  PRINT  'THIS  PROGRAM  CALCULATES  THE:  HOURLY  TRANSIENT* 

A  PRINT  ’HEAT  FLOW  THRU  A  MULTILAYERED  RUH  R J  NG  COMPONENT* 

5  PRINT  'EY  USING  THE  < TRANSFER  FUNCTION  METHODS  * 

6  PRINT  * PRESENTED  BY  ASHRAE . HANDBOOK  OF  FUNDAMENT ALS  ( 1 v?V ) , “ 

7  PRINT 

8  PRINT 

9  PRINT  'YOU  WILL  BE  PROMPTED  TO  SUPPLY  NEEDED  INFORMATION.1 

TO  PRINT 

11  PRINT 

15  PRINT “THIS  PROGRAM  WILL  PRINT  OUT  A  ITERATIONS  OF  HEAT  FLOW  FOR* 

16  PRINT 'EACH  SET  OF  DA 'A  ENTERED.  IF  YOU  WISH  MORE  OR  LESS  ITERATIONS* 
:I.7F‘RINT  “  CHANGE  STATEMENT  NUMBER  410.  * 

18PRINT 

1 9PRINT "  ***************** ###*###****#**** Y*****##***##*###*##*##  1 

20PRINT 

21  PRINT 

26  DIM  T  < 1 00 ) »  A ( 1 00 ) *  B ( 20 ) >  D i 20 ) *0(100) .01 ( 100) . T 1 ( 20  > 

28PRINT  *  HOW  MANY  HOURLY  SURFACE  TEMPERATURES  DO  YOU  HAVE ’ ; 

30  INPUT  I 

31  PRINT 

32  PRINT 

40  PRINT  ‘HOW  MANY  B-COEFFICIENTS  DO  YOU  HAUF.  “  i 
50  INPUT  II 


51 

F,0 


pm 


60  PRINT  '  HGU  MANY  It-COE.FF IC.TENTS  DO  YOU  HAVE  *  i 
70  INPUT  12 
71  PRINT 
72  PRINT 

80  PRINT  'WHAT  IS  THE  VALUE  OF  YOUR  C-COEFFICIENT  *  t 

90  INPUT  C 

91  PRINT 

92  PRINT 

100  PRINT  ’WHAT  IS  YOUR  AVERAGE  INDOOR  TEMPERATURE  IN  DEGREES  FAHRENHEIT* 

110  INPUT  T2 

111  PRINT 

112  PRINT 

115  PRINT  "TYPE  IN  EACH  HOURLY  TEMPERATURE’; 

120  MAT  INPUT  A 
122  PRINT 

125  PRINT  ’  TYPE  IN  EACH  B-COEFFICIE'NT  *  f 

126  MAT  INPUT  B 

129  PRINT 

130  PRINT  *  TYPE  IN  EACH  D-COEFF ICJENT * t 

131  MAT  INPUT  D 

140  FOR  N=I 1  TO  < Ilrl-l  ) 

150  S=N-I 1+1 
160  T  <  N>  =A ( S ) 

170  NEXT  N 

171  PRINT 

172  PRINT 


It 


I8*>  t-  h  ,  M  •  *  THE  EOI.  TUNING  ARE  THE  bURFAGE  I  EMPERATUKFS  U:>Hi  l  N  )  H T !•;  PRObEnri 
I  H  i  PRINT 

tft^-F'fclNT  "  HOUR  TRIP* 

183  PRINT 

190  FOR  N~ 1 1  TO  <  T  t  f  T - 1 ) 

191  SI.  =N-1 1  +  .1 
200  PRINT  SlfT(N) 

?10  NEXT  N 

220  FOR  H=  1  TO  <  X  :l  - 1 ) 

230  L  ~ii+  T 
240  T(M)=T<L) 

250  NEXT  Pi 
27 7  PR  INC 

27 a  PRINT “THE  B-COEFF TCI ENTS  ARE* 

279  PRINT 

280  FOR  J= 1  TO  11 

301  PRINT  USING  307, J,R( J) 

302:  t*t  tt.ltftH 

303  NEXT  J 

332  FOR  Ul  =  l  TO  13 

333  M1=U1-1 

334  TKHl  )=D(IJ1  ) 

330  NEXT  U1 

337  PRINT 

338  PRINT  ’THE  D-COEKFICIENTS  ARE* 

339  PRINT 

340  FOR  K-0  TO  (12-1) 

36  1  PRINT  USING  362pK»D<K> 

362 J  ***  #*. ****** 

370  NEXT  K 

371  PRINT 

372  PRINT*  HOUR  HEATFLOW  BTU/HR  * 

373  PRINT 

410  FOR  A*1  TO  4 

420  FOR  P=I.l  TO  <1  +  11-1  ) 

430  X=P-<I1-1) 

450  FOR  S=1  TO  II 
460  P1®<P+1)-S 
470  T1(S)~T(P1) 

480  NEXT  S 
490  B2=0 

500  FOR  U“1  TO  II 
510  B=B(U)#T1 <U) 

520  B2=B2+B 
530  NEXT  U 
540  G  <  0 ) =0 
550  D2=0 

560  FOR  V±1  TO  <12-1 > 

570  01=P-0 

600  B=D<O)*CJl<01 ) 

6.10  D2=D2  +  D 

620  NEXT  0 


